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It is shown that, by imposing reparametrization invariance, one may derive a variety of stochastic 
equations describing the dynamics of surface growth and identify the physical processes responsible 
for the various terms. This approach provides a particularly transparent way to obtain continuum 
growth equations for interfaces. It is straightforward to derive equations which describe the coarse 
grained evolution of discrete lattice models and analyze their small gradient expansion. In this 
way, the authors identify the basic mechanisms which lead to the most commonly used growth 
equations. The advantages of this formulation of growth processes is that it allows one to go 
beyond the frequently used no-overhang approximation. The reparametrization invariant form 
also displays explicitly the conservation laws for the specific process and all the symmetries with 
respect to space-time transformations which are usually lost in the small gradient expansion. 
Finally, it is observed, that the knowledge of the full equation of motion, beyond the lowest order 
gradient expansion, might be relevant in problems where the usual perturbative renormalization 
methods fail. 
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I. INTRODUCTION 



There has been a great deal of recent work on the formation, growth and geometry of interfaces (Family and 
Vicsek, 1991; Mcakin, 1993; Barabasi and Stanley 1995; Halpin-Healy and Zhang, 1995). These studies are relevant 
to a variety of experimental situations including biological growth (Eden, 1958), the propagation of flame fronts 
(Sivashinsky, 1977, 1979), fluid flow in porous media (Cieplak and Robbins, 1988; Martys et al., 1991) and atomic 
deposition processes (Tu and Harris, 1991; Messier and Yehoda, 1985; Bales et al. 1990; Tang et al. 1990) such as 
the technologically important molecular beam epitaxy (MBE). On a more fundamental level, some of these processes 
are prototype of far-from-equilibrium physics without a Hamiltonian formulation (Hohenberg and Halperin, 1977). 
Recent advances have shown that it is nevertheless possible and useful to categorize the systems into universality 
classes. There have been two principal approaches for the theoretical analysis of such problems. The first (see, e.g., 
Meakin, 1993) is based on computer simulations of discrete models and often it provides useful links between analytic 
theory and experiments. The second approach (Krug and Spohn, 1990; Kardar, 1994; Halpin-Hcaly and Zhang, 1995) 
is to describe the dynamical process by stochastic differential equations. This procedure neglects the short length-scale 
details but provides a coarse grained description of the interface which is suitable for characterizing the asymptotic 
scaling behavior. There are two essential steps to be carried out: first, one must deduce the continuum equation, 

2 



and second, in order to obtain the scaling behavior the equation needs to be solved or analyzed by renormalization 
group techniques (Ma and Mazenko, 1975; Forster et al., 1977; Wolf, 1991; Frei and Tauber 1994; Sun and Plischke, 
1994). The principal theme of this paper is the derivation of the stochastic partial differential equation appropriate 
to the different physical processes responsible for the growth. We will show that the principle of reparametrization 
invariance (R-Invariance) can be used in a straightforward, yet powerful, manner to derive the continuum equation. 
This approach may also be useful in other contexts such as the study of avalanche dynamics associated with an 
interface moving in a random medium (Makse, 1995) and for the study of interfaces in complex dynamical systems 
(Kapral et al, 1994). 

We restrict ourselves to local growth processes in which the growth rate is a function of the local properties of 
the interface. (Non-local effects are discussed in the last section in terms of partial differential equations for two 
coupled fields.) Traditionally, continuum equations have been derived in the no-overhang approximation (Mongc 
representation). The height of the interface on a reference substrate plane, h(f,t), is assumed to be a single valued 
function of the lateral coordinates r and time t. While this assumption is valid and yields the correct scaling behavior 
in many instances, it has been demonstrated that some essential details are left out in this approach (Keblinski et 
aZ.,1994, 1995, 1996). For example, in the columnar growth regime of sputter deposition of thin films, the merging and 
regeneration of columns arises due to the presence of outward flaring columns. This phenomenon, observed even with 
a normal incidence of the depositing atoms, is not captured within the no-overhang approximation. The R-invariance 
principle will allow us to derive the continuum growth equations without such restrictions - the no-overhang situation 
can then be obtained as a special case. 

Previously, continuum growth equations have been derived either directly from discrete models by Vvedensky et al., 
(1993) or by using methods based on preserving the symmetries and conservation laws of the system (see, e.g., Hwa 
and Kardar, 1989; Kardar, 1994). The latter approach has been widely used to identify universality classes of discrete 
atomistic models (Liu and Plischke, 1988; Huse et al., 1990; Racz et al., 1991; Das Sarma and Tamborenea, 1991; Das 
Sarma and Ghaisas, 1992; Kotrla et al., 1992; Amar and Family, 1993; Krug et al., 1993). While this method often 
unambiguously yield the lowest order terms in a gradient expansion that could be present in the stochastic differential 
equation, physical considerations have to be invoked to decide whether such terms are indeed present or not. In the 



next section we discuss previous work on the derivation of the dynamical equations. Section [II describes how the 
R-invariance principle can be used to obtain the growth equations in a variety of physical situations. The last section 
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summarizes our results. 

The scope of our article is limited the emphasis is on theory and exclusively on the derivation of continuum growth 
equations and discussions regarding their geometrical and physical content. Here and there, we shall mention results 
which have been obtained elsewhere for the solutions of some of these equations and their relationship to physical 
processes such as MBE or discrete models. Our discussion will certainly fail to be complete, a task which goes far 
beyond our goal of providing a self-contained exposition of the derivation of growth equations. We apologize to the 
many authors whose work is not appropriately mentioned in this article. 

A quite complete survey of this field can be found in several review articles, each with its own focus. Among 
the early reviews, the book by Family and Vicsek (1991) and the contributions by Meakin (1988a) and Krug and 
Spohn (1990) have become standard references. Specifically, the first contains a collection of important reprints with 
a commentary, the second focuses mostly on the phenomenology of simple models of growth while the third contains 
a discussion of analytical approaches. Another early review is the book by Pietronero and Tosatti (1985) which 
describes fractal growth processes. The effects of a random substrate on wetting phenomena have been considered 
among others by Pfcifcr et ai, (1989) and Giugliarelli and Stella (1991, 1994). Recent work by Halpin-Hcaly and 
Zhang (1995) concentrates on the Kardar-Parisi-Zhang (1986) equation and reviews the efforts of the last decade 
devoted to elucidating its scaling behavior. From this central theme, it extends to cover the field in a pedagogical 
manner, with special emphasis on directed polymers in random media. Numerical simulations of discrete growth 
models are nicely reviewed in the work by Meakin (1993). A quite complete elementary book with an extensive list of 
references by Barabasi and Stanley (1995) has recently been published, where a lucid discussion of interface growth 
in random media, including a summary of recent developments, can be found. Finally, a detailed discussion of the 
phenomenology and some theoretical approaches for the study of MBE growth can be found in the book by Tu and 
Harris (1991) (see also Tang and Nattermann, 1991; Villain et ai, 1992; Kcsslcr et ai, 1992; Hunt et ai, 1994; Pal 
and Landau, 1994; Siegert and Plischke, 1994; Das Sarma, 1994; Das Sarma et at, 1996). 

II. PREVIOUS APPROACHES 

In this section we will briefly review the traditional method for deriving local growth equations. With very few 
exceptions (Meakin et ai, 1986; Maritan et ai, 1992), this is accomplished within the familiar no-overhang approx- 
imation (the Monge representation). Denoting by h(x,t) the single-valued interface height function of the lateral 
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coordinates x and time t, one may generally write a growth equation in the form 

WM = G[h%t)]+rii2,t) (1) 

where G is the deterministic growth term and rj is the noise. This is essentially Newton's law of motion in which 
the inertial term has been neglected in comparison with the dissipative force. The neglected term is irrelevant in 
determining the scaling behavior in the asymptotic regime. This regime, in Fourier space, is related to the behavior of 
correlations functions in the limit of small frequencies cu — > 0. While the left hand side (LHS) of Eq. ([!]) is proportional 
to u>, the inertial term is of order oj 2 and hence it is irrelevant. 

Consider a D dimensional substrate of size L (the physical case corresponds to a L x L substrate with D = 2) 
and then define the mean height of the growing film and its roughness by 

h(t,L) = Jb I d D xh(x,t), (2a) 



1/2 

W(t,L) = (Jd J d ° x [HS,t)~h(t,L)] 2 ^ (2b) 



where (. . .) denotes an average over different realizations of the noise (samples). 

Starting from a flat interface (one of the possible initial conditions) it was conjectured by Family and Vicsek (1985) 
(see also Plischke and Racz ,1985, Family and Vicsek, 1991) that a scaling of space by a factor £ and of time by a factor 
l z (dynamic scaling), rescales the roughness W by a factor £ a , as appropriate for a self-affine surface (Mandelbrot, 
1986), i.e., 

W(£ z t,£L)=£ a W(t,L), (3) 

which implies that 

W(t,L)=L a f(t/L*). (4) 

If, for large t and fixed large L (t/L z — » oo), W saturates, then f(x) — > constant as x — > oo. However for fixed large 
L and 1 <C t <C L z , one expects that correlations of the height fluctuations are set up only within a distance t x / z and 
thus W must be independent of L. This implies that for 

f(x) ~ x 13 with P = a/z. (5) 

Thus dynamic scaling postulates that: 



W(t, L) ~ t 13 1«(«L 2 , 



(6a) 



~ L a i>L z . (6b) 

The roughness exponents a and the dynamic exponent z characterize the self-afhne geometry of the surface and its 
dynamics respectively. 

The above considerations do not include short time behavior. In that limit one may expect that the height 
fluctuations are uncorrelated and therefore h(x,t) — h(t) behaves as f T](T)d,T where rj(r) is a white noise with zero 
average and (?7(t)?7(t')) cx <5(t — r'). This leads to W(t, L) ~ i 1 / 2 independent of L and is the random deposition 
regime. Note that this dynamic scaling is not related to the standard dynamic scaling in critical phenomena. The 
latter involves correlations which are space time separated whereas here they concern the behavior at the same time. 
Furthermore, after the interfacial roughness has saturated (t 3> L z ), one may expect that correlations of fluctuations 
6h(x, t) = h(x, t) — h(t, L) such as: 

C L (x,t;x* ,t') = ([5h(x,t)-Sh(x",t')} 2 ^ (7) 
show the more traditional dynamic scaling: 

C L {x, t + r; a* t) = \x- x'\ 2a g L \ nz , , ] -^±) (8) 



L 

with z' not necessarily equal to z. Note that the exponent a must be the same as before since 

1 



d u xd u x'C L {x, t; x', t) = W 2 (t, L) ~ L 2a (t » L z ). (9) 



L 2D 

The key idea behind the traditional derivation of G and r\ in Eq. (|1|) is the identification of the symmetries and 
conservation laws of the system (Kardar, 1994). The deterministic term G is expanded in powers and combinations 
of the h field and the relevant (in the renormalization group sense) lowest order terms that are consistent with the 
symmetries and conservation laws are retained. The lowest order contribution is a constant Go which can always be 
set to zero by the transformation h — > h + G$t. This amounts to the use of a reference frame that co-moves with 
the interface. Symmetries that need to be considered include invariance under translation in time and space, and 
along the growth direction, rotational and inversion symmetry about the growth direction. For example, translational 
invariance in the growth direction, which holds in most practical cases, rules out terms proportional to powers of h. 
An explicit dependence of G on time or position x contradicts time or space translational invariance, which usually 
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holds, and can therefore be neglected. In addition it is helpful to classify terms which might occur in the expansion 
in two categories: those which conserve the number of particles and those which do not. To derive the equation of 
growth for a conservative dynamics, the key observation (Villain, 1991; Kardar, 1994) is that the deterministic part 
must have the form of a continuity equation: 

^ (10) 

where the macroscopic current j(x,t) describes the flux of atoms on the surface. The current j(x,t) arises in general 
from differences in the local chemical potential /i(x,t). The relation between j(x,t) and n(x,t): 

j(x,t) = -V f i(x,t) (11) 

describes the fact that atoms drift to regions of minimum chemical potential. The simplest source of chemical potential, 
which is gravitational energy \x oc h, led Edwards and Wilkinson (1982) to formulate the equation 

^jM = uV 2 h(x,t)+ V (x,t). (12) 

The noise here arises due to the stochastic nature of the incoming flux of atoms and, as such, it is non-conservative. The 
noise can also result from the probabilistic nature of the diffusion process on the surface or from thermal fluctuations. 
In the latter cases, the number of atoms on the surface remains the same and the noise will be called conservative. 
The noise correlation functions are different in the two cases and are given by: 

{r](x, i)j7(f , t 1 )) = 2D 6 D (x - tf)6(t - t') non-conservative (13) 
(r){x,t)r ) (x',t')) = {-DiV 2 + D 2 V*)5 D {x-x')5(t-t') conservative (14) 

We shall not consider the case where the statistics of -q is not gaussian (see e.g. Zhang 1990, and Krug 1991, Horvath 
et al, 1991, Lam and Sander 1993) or that in which it is long range correlated in space or in time (see Medina et al. 
1989, Peng et al. 1991, Lam et al. 1992). 

Usually, the effect of gravity on deposition processes is totally negligible with respect to other effects. However, 
the Edwards- Wilkinson equation also applies to other situations and mainly to the case where atoms are allowed to 
evaporate from the surface. This effect, which does not conserve the number of particles, is nevertheless described 
by a term v\/ 2 h in the growth equation. We repeat here the derivation, given by Villain (1991). In general, if 
the deterministic growth term G describes dynamics which minimizes a potential V[h] (which is a functional of the 
interface configuration h(x,t)) it will have the form 



6h(x,t) ^ 
where 5/Sh(x,t) denotes a functional derivative. The effect of evaporation is to minimize the surface area. Since the 
excess surface area due to roughness is given, in a small gradient expansion, by V[h] = J d d x (v/i) (Br uinsma and 
Aeppli, 1984), Eq. ( [l5| ) directly gives the first term on the right hand side of Eq. (|l2|). It is interesting to note already 
at this stage how a relevant feature of the process, i.e., whether the number of particles is conserved or not, is lost 
in the small gradient expansion of the growth equation. A derivation of the growth equation from the R-invariance 
principle avoids such ambiguities. 

Let us return to conservative equations and look for the simplest equation one gets by neglecting gravity. ^From the 
above discussion, it is natural to expect that the next term in the expansion of the chemical potential will be a term 
proportional to the local curvature fj, cx V 2 /i. This term, via Eqs. (|l0|), (|TT|), will favour again a flux of particles away 
from the local maxima towards the local minima. Hereafter, the resulting relaxation mechanism will be called surface 
diffusion (Villain 1991; Siegert and Plischke, 1992, 1993; Das Sarma et ai, 1996). Since this dynamics constrains the 
atoms to stay on the surface, this relaxation will be much slower than the one provided by evaporation. The resulting 
linear equation is 

?MM = - K V*h(x 1 t)+ v (x,t) (16) 
ot 

and the same considerations discussed previously apply to the noise term here. The two equations can be combined 
into one which accounts both for gravity and for surface diffusion. In experimental situations of MBE, however, the 
typical coefficients v and k are such that the effects of gravity are relevant only on length scales much bigger than 
the typical sample size. We note that both terms discussed so far satisfy all the symmetry properties mentioned 
previously. Note that while rotational invariance in the substrate plane is satisfied by these terms, the full space 
rotational invariance is lost. This is a consequence of the small gradient expansion in h rather than a physical 
characteristic of the processes. For example, atoms diffuse on the interface in the same manner irrespective of the 
direction we assign to the z axis. Again, we will recover the full space rotational invariance, and distinguish terms 
which satisfy this symmetry and those that do not, once the R-invariant form of the equation is derived. 

The lowest order non-linear term which one could include in G leads to the Kardar-Parisi-Zhang (KPZ) equation 
(Kardar et al, 1986) 

= vV * h ^ t ) + ±(Vhf + V (x, t). (17) 



Here the mechanism of relaxation of surface fluctuations is the same as that of the EW equation The origin 

of the non-linear term lies in the driving force of the deposition process that is perpendicular to the interface (or 
lateral growth in the discrete ballistic deposition model - for a review of ballistic deposition, see Meakin, 1993. 
See also Ko and Seno, 1994 for an off-lattice simulation). While the simplest growth term, a constant Go, can be 
eliminated by choosing a co- moving frame on the interface (sec above), the non linear term of the KPZ equation 
cannot. The resulting dynamical process is intrinsically irreversible. Interesting generalizations of the KPZ equation 
to a multi-component model have been studied by Doherty et al. (1994). 

Among the notable properties of the KPZ equation are its relations with a large number of other problems, including 
the statistics of directed polymers in random media (Kardar and Zhang, 1987, Fisher and Huse, 1991, Kim et al., 
1991, Halpin-Healy and Zhang, 1995) and the Burgers equation for fluid dynamics (Forster et al., 1977). The Galilean 
invariance in the latter problem has profound consequences on the properties of the KPZ equation, and leads to the 
exponent identity z + a = 2 (Forster et al., 1977; Meakin et al., 1986). This invariance translates the problem of 
interface growth into an invariance of the equation for an infinitesimal tilt of the substrate plane. This invariance can 
be seen as a remnant of the full D + 1 rotational invariance of a growth process occurring in the direction normal 
to the interface. This feature can be fully appreciated in the R-invariant form (Maritan et al., 1992) of the equation 
which describes an isotropic growth process driven by a pressure or occurring from the condensation of a vapor. In 
the gradient expansion leading to Eq. (|i"7j), rotational invariance is retained only for infinitesimal transformations. 
However, as we shall see, there are other mechanisms, which are not rotationally invariant, which lead to a growth 
equation of the form in the small gradient expansion. 

Note that in the KPZ equation the h — > —h symmetry is broken: there is a definite growth direction. Also the KPZ 
equation describes processes in which the number of particles is not conserved. The (Vft.) 2 term cannot appear in a 
situation, such as Molecular Beam Epitaxy (MBE), where conservation of the number of particles on the interface is 
expected to hold. 

One may invoke dynamic scaling to show that, in general, an equation of the type: 

°-g = (-l)^K n (V 2 ) n h + r 1 with (r 1 (x,t) V (x>,t'))=2Dr m 6 D (x-x>)6(t-t>), 

where m takes on the two values or 1 with To = 1 and Ti = —V 2 , has the exponents 

2(n-m)-D 

z — 2n and a = . 

2 
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Thus for the Edwards- Wilkinson model Eq. ( |l2| ) and for Eq. (|l6|), one obtains (to = in both cases) 



2-D 



and 



A-D 



z = 2 



a = 



z 



= 4, 



a = 



2 



2 



respectively. The situation with a < corresponds to W(t — > oo,L) ~ constant (a = leads to logarithmic 
corrections) . Since the continuum equation was derived in a power series of gradients of h with only the linear term 
being retained, the result a > 1 is a matter of concern, since local slopes on a distance of order i scale as \Vh\ ~ l a ~ l . 
This suggests that neglected non-linear terms need to be considered as well (for a discussion on self-similar (a = 1) 
and self-affine surfaces see Mandelbrot, 1986 and Meakin et al, 1986). The requirements that the basic symmetries 
be obeyed and that in the surface diffusion dominated regime the continuity equation ( |To| ) be satisfied lead to: 



in which all the terms up to fourth order in the gradient expansion of h have been retained. Physically v is set to 
zero since it arises from a gravity-like chemical potential which, as mentioned above and as we will discuss below, 
is negligible. The requirement that the surface current j results from a chemical potential would rule out the A2 
term, even though this term is more relevant, from a renormalization group point of view, than the Ai term. It will 
be clear, from the R-invariant form of the equation, that the Ai term arises as the second one in an expansion of a 
gravity or surface tension term, the first one being the v term. It is interesting to note that recently (Das Sarma and 
Kotlyar, 1994; Kim and Das Sarma, 1995) it has been proved that, even if the coefficient of the Laplacian is zero 
in the starting equation, a v > is generated under renormalization by the A2 term, leading ultimately to the EW 
behavior for surface fluctuations. 

Recently a growth equation has been derived by Vvedensky et al, (1993) starting from a microscopic solid-on-solid 
model with deposition, desorption and diffusion through a master equation for the surface dynamics. In the continuum 
limit, the growth equation is obtained - desorption leads to non-zero v and A values while surface diffusion produces 
the k term. 

Similar considerations have been used to derive growth equations for interfaces in a disordered medium (Koplik and 
Levine, 1988; Jensen and Procaccia, 1991; Parisi, 1992; Nattermann et ai, 1992; Jiang and Hentschel, 1992; Sneppen, 
1992; Nolle et ai, 1993; Olami et ai, 1994; Galluccio and Zhang, 1995) such as domain walls in random field Ising 
models (Nattermann, 1985; Fisher, 1986; Ji and Robbins, 1992; Narayan and Fisher, 1993) and the invasion of one 
fluid into another within a porous medium (Nittmann et al, 1985; Stokes et al., 1988; Buldyrev et at, 1992; He et 



dh(x, t) 
dt 



zA7 2 /i(x, t) - n\7 4 h(x, t) + AiV 2 (V/i) 2 + \2V[Vh(Vh) 2 } + n(x, t), 



(18) 
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at, 1992; Tang and Leschhorn, 1992; 1993; Leschhorn, 1993; Delker et at, 1996). The disorder serves to pin parts of 
the interface. This feature is captured by replacing the thermal noise with a quenched random noise generated by the 
disorder rj(x, h) governed by correlations of the form 

(t](x ) h)ri(^,h')) =2D a 5 D {x-x')A(h-ti) (19) 

where (. . .) represents an average over different realization of the randomness and the function A characterizes the 
nature of the quenched disorder. It is important to note that disorder will break the translational invariance in the h 
direction. Hence a constant Go term can no longer be eliminated by choosing the co-moving frame, since under the 
transformation h(x, t) — > h(x, t) + G$t, it will reappear in the argument of the noise r)(x, h) — > r)(x, h + Got). Here, 
Go plays the role of the driving force: for small Go the interface will find sooner or later a surface where the values of 
the pinning forces are strong enough to inhibit further growth. In contrast, if Go is very large, h in the argument of 
the noise in the co-moving frame can be neglected with respect to Got, and one expects to recover the behavior of the 
corresponding dynamics without disorder, i.e., A(h — h') oc S(t — t'). In practice the interface will move so fast that it 
will sample so many values of the disorder in a small time interval, that the overall effect is that of a time dependent 
noise with correlations of the form (|l3|). We shall not deal in the following with surface growth in disordered media, 
a steadily growing subject which is reviewed in the book by Barabasi and Stanley (1995). 

III. EQUATIONS FOR A GROWING SURFACE AND REPARAMETRIZATION INVARIANCE 

The most general Langevin equation for the evolution of a surface in a D + 1 dimensional space has the form 

d t r{s,t) = n{s,t)g[f{s,t)]+F{s,t), (20) 

where the D + l dimensional vector f(s, t) = {r Q (s, t)}^^ 1 runs over the surface as s = {s l }^ =1 , varies in a parameter 
space. (See Appendix ^ for a brief summary of the elements of differential geometry). In Eq. (^0|) n stands for the 
versor normal to the surface at r while Q contains a deterministic growth mechanism that causes growth along the 
normal h to the surface and is a functional of r itself. F is a random force acting on the surface. Eq. (^0|) derives 
from Newton laws in the limit of a massless surface when the inertial term Off can be ignored with respect to the 
dissipative force. Note that the time derivative of r has to be parallel to the normal to the surface. This is because 
t = s° can be regarded as the D + 1 th coordinate, and (s, s°) is a curvilinear coordinate system. If the growing surface 
invades the D + l dimensional space, this parametrization is legitimate since the metric tensor is positive definite. 
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However s° is the absolute time, and changes of parametrization cannot involve this variable. This is satisfied only if 
the elements go,i = dtf ■ dif vanish, which implies that dtf _L dif for all i = 1, . . . , D, and therefore dtf || n. 

Independent of the physical mechanisms entering the various terms of Eq. (p0|), which specify the form of Q and 
the properties of F, this equation has to satisfy the fundamental requirement of reparametrization invariance (R- 
invariance) . This requires that only quantities that are independent of the choice of the parametrization s, such as 
those referring to the local geometry of the surface, like the curvature, can appear in the equation. As with any other 
symmetry, reparametrization invariance poses constraints on the possible forms that Eq. (|2^) can take. 

The basic elements of differential geometry used in the derivations and some calculational details are presented in 
the appendix. 

It will be convenient to express the deterministic part of Eq. (|2(i|), as well as the noise, as a sum of different terms 

G = G a +Qb + ■■■ and F = F a +F b + .... 

The derivation of the equation can then be split into those of the individual terms that are expected to appear in a 
given physical system. 

It is often convenient and sufficient to describe the interface in the Monge form, i.e., with f — (x, h(x)). Henceforth, 
at variance with section [n], we will use the notation x instead of x to stress that the Monge form is one of the many 
parametrizations one can use. h(x) is the coordinate in the direction normal to the substrate and since it must be a 
single valued function of x the interface cannot have overhangs. In this form ( pp| ) becomes 

d t h(x,t) = ^g(g + V ) (21) 

where g = 1 + (V/i) 2 is the determinant of the metric tensor and the relation n = (— V/i, l)/y / 5 has been used (see 
Appendix ^|). Q and rj = F ■ n, as before, are the amplitude of the deterministic force and the noise in the normal 



direction respectively (note that dth(x,t)/y/g is the normal velocity n ■ dtf - this is easy to prove using Eq. (A14)). 



The discussion of the noise term is presented in Section [II. B 



A. Deterministic evolution 

As stressed before, reparametrization invariance requires that Q depends only on intrinsic geometric properties of 
the interface such as the mean curvature H or, when Q is not rotationally invariant, on scalar products of the normal 
fi with some fixed vector v. 
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The physical meaning of Q is particularly evident in the case in which it can be derived from a potential: 

_ , . . i sniffs)] 

where n is an R-invariant functional of f. The 1/Jg term in Eq. (|22|) appears since 1 5 /}r~\ — S ^~— I is the 
R-invariant Dirac's delta function (see Appendix A). In this case, the dynamics tends to minimize the potential energy 



n of the surface. Moreover, if the random force is properly chosen (this will be discussed in section [II.B), i.e., if it 
is not conservative, the system approaches a steady state whose distribution of f is given by exp{— [3n[f}} where (3 is 
related to the correlations of the noise (Hohenberg and Halperin, 1977). 

The R-invariance of n guarantees that the functional derivative in Eq.(|2^) is a vector parallel to the normal as 
required by Eq. (^0|), and the 1/-Jg factor guarantees R-invariance of the functional derivative. Indeed R-invariance 
of n implies that n[f(s')] = n[f(s)] for any reparametrization s'(s). For an infinitesimal transformation s' = s + e(s), 

P(s) = f(s') = f{s) + e(s) ■ df(s) + 0(e 2 ) 

so that 

n[f{s!)] = n[f(s)] + [ d D se i {s)dif{s) ,)H 



Sf(s) 

Since dif is a vector in the tangent plane of the surface, the second term in this equation vanishes for any function e 
only if the functional derivative is parallel to the normal n . 

In the Monge form, if Q derives from a potential n we find (see appendix |b|): 

G = - n -vli>f = -ih- (23) 

We now proceed to consider the simplest possible terms that can appear in Eq. (|2^). We first give the expression 
in a general parametrization and then discuss the Monge form and the expansion in small gradients of h(x,t). 

1. Surface tension 

The simplest physically motivated term in the Hamiltonian of a surface is proportional to the total area A — 
J d D Syfg and produces a force that tends to minimize the surface area. This term is usually referred to as arising 
from surface tension. The functional derivative yields (see Appendix ^) 



(24) 
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where A is the Beltrami-Laplace operator defined in Eq. (A4). For 7i s = v s A, using (A7), one gets: 



Q s = v s h ■ Ar = v s H, 



(25) 



where H is the mean curvature. 



Using Eqs. (p3|), ( |25| ) and 7i s = v s § dPx^f\ + (V/i) 2 , one finds, in a small gradient expansion in the Monge 
parametrization (Eq. (|2l|)): 



W 2 h(x,t) + ^V 2 /i(V/i) 2 - iv • [V/i(V/i) 2 ] 



(26) 



The main physical mechanism which produces a term like this is evaporation. Following Villain (1991) we observe 
that the evaporation rate will be proportional to the difference between the chemical potential of the solid fj, s and 
that of the vapor fj, v . On the surface, the former will depend on the local geometry and hence: 

<9 t r| eva p. = nB(n s (r) - fx v ). 



We can then expand fi s in powers of the local curvature. Along with the zeroth order term, corresponding to the 
chemical potential for a flat surface, there will appear a term proportional to the curvature H as in the r.h.s of Eq. 
25|>. 



2. Pressure 

The second simplest geometrical property on which Ti may depend is the volume enclosed by the surface, which 
may be written as V = 55W / d D sy/g(s) r (s) • n(s). A linear dependence of Ti. on the volume physically represents a 
pressure term. If Ti v = —XV, the pressure A > encourages an increase in the volume, while if A < the force in Eq. 
( p2| ) acts to deflate the surface. The infinitesimal volume variation on the surface element da = d D s^Jg is given by 
dan ■ Sr, so the functional derivative of TL V in Eq. (|23) gives: 



Gv = -n ■ — —^7- = A, 
V9 or 



(27) 



that in the Monge form becomes 



VdSv = AV5 = * 



l+2^) 2 



(28) 



Equation (p0|) with Q = Q s + Q v is one of the R-invariant forms (Maritan et at, 1992) of the KPZ equation: 
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d t r(s,t)\ det =n(v s H + X). (29) 

which in the Monge representation reads: 

V/i 

d t h(x,t)\ det = Vay/gS. ' == + X^9- (30) 

Indeed, to lowest order in the gradient expansion, Eq. (|l7]) is recovered (apart from the constant term A which can 
be absorbed by redefining h — ► h + Xt). The complete R— invariant KPZ equation derives from the Hamiltonian: 

H KPZ = /" dPx(i> ay /g-\h), (31) 

which, however, is unbounded as h — > oo. Thus, even with a suitable noise term, Eq. (|3^) does not have exp(—07i, KPZ ) 
as the equilibrium distribution of h. This is not surprising and is related to the presence of a pressure that makes the 
system grow forever. In a sense, the interface growth is intrinsically irreversible. Note that A couples only to the k = 
mode of ft. in a Fourier expansion. This would suggest that steady state correlation functions, such as ([h(x) — h(y)] 2 ), 
are independent of A. A derivation of the KPZ equation from the functional derivative of a free energy with a volume 
and a surface term, was also obtained by Grossmann et al. (1991) in a more complex way. An alternative derivation 
was also given in Kcblinski et at, (1996). 

3. Curvature energy 

The potential TL may also depend on the curvature H of the interface. In general, this dependence can be expressed 



in a power series expansion (the zeroth order term has already been considered in Section III.A.l) 



H c = J d D s^(K 1 H + k 2 H 2 + ...) = H c ,x + H c ,2 + ..■■ (32) 

The physics behind the first term reflects the difference in the mechanical properties of the media divided by the 
interface. Indeed for k± > 0, large negative curvatures are encouraged while positive ones are depressed. The 
functional derivative is carried out in the appendix |D| with the result 

e C)1 = -^n.^ = K J j IF - > AH ■ (•»•'-!) 



y/9 Sf 

where Xi are the eigenvalues of the matrix of the coefficients of the second fundamental form and express the principal 
curvatures of the surface. Since H = X\ in D = 1, Q Ct i vanishes. This is a consequence of the Gauss-Bonnet theorem 
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that states that the integral of the Gaussian curvature K on a closed surface is a constant. Since H = K in D = 1, 
the variation of TL c ,i is zero. 

In terms of h(x), Eq. (|33|) takes the form: 



(34) 



In the small gradient expansion no linear term in h arises. The curvature term, as expected, breaks the symmetry 
h — > —h in the growth equation. 

Similarly, higher powers of H are easily worked out ( see Appendix |d|) . For the p th term in Eq. (^) we find: 



i . sn t 

n ■ — — 

V9 ° r 



(35) 



This equation generalizes Eq. (25), that is the p = case, and Eq. ( |33| ) (p = 1). The effect of these terms is more 
transparent in the Monge parametrization. Of particular interest is the K2 term in the expansion of the curvature 
potential: 



Vg5c,2 = -2n 2 (y 2 ) 2 h- 



(36) 



4. Rotational Invariance 



All the C/'s we have considered so far are invariant under rotations in the D + 1 space. This implies (unless the 



noise term breaks this invariance) that the evolution equation in its form Eq. (21) is invariant under the following 
(infinitesimal) symmetry transformations: 

x' = x— eh(x, t), 



and 



h'(x',t) = h(x,t) + e-x. 

Indeed, it is straightforward to verify that the normal velocities dt^fe*) an( j &th (x ,t) equal. 

ygfet) y/g'(x't) 

In the small gradient expansion of Eq. ([50|), Eq. (|l7|) is obtained after redefining h h + Xt. As a consequence, the 
above symmetry transformations become 

x' = x — eXt, 
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and 

h'{x',t) = h(x,t) +e-x. 

where —eh(x, t) has been omitted in the first equation, since it contributes only to the higher order powers neglected 
in Eq. (|30|). These symmetry transformations are exact for the KPZ equation and can also be deduced from the 
Galilean invariance of the related Burgers' equation (Forster et ai, 1977; Halpin-Healy and Zhang, 1995). 

5. External potentials 

In the presence of an external potential, 7i depends explicitly on r. The simplest example is that of a gravitational 
field. The variation 8TL g of the gravitational energy, when f — > r + 5r, is proportional to the variation of the mass 
pdP s^fgn ■ 5r {p is the mass density) multiplied by the acceleration a g of the gravitational field and the "height" z ■ r, 
where z is the direction of the gravitational field. Then 

5n g = v 9 [ d D Syfg{f- z)h ■ Sr. (37) 



where v g = pa g . It is easy to see that such a term breaks translational invariance, i.e., the growth equation changes 
if r(s) — > r(s) + Tq. This generally applies to any potential V(r) for which Ti = j d D Sy/gV(r), unless V(r) = v s 



which yields a surface tension term (see [II. A. 1), or V(r) = —Xn ■ r/(D + 1), corresponding to the case of a pressure 



potential dealt with in [II. A. 2. Since translational invariance is expected to hold, these terms are not considered. In 
the Monge parametrization, the gravitational energy, proportional to J d D xh 2 (x) , leads to a linear term in h in the 
equation of motion. 

6. Orientational energy 

The simplest, translational invariant term which breaks rotational invariance, results from considering a potential 
which depends on the local orientation of the surface: 

n z = - J d D s^ X (n z ), (38) 

where n z = h ■ z, z is some fixed direction in the D + 1 dimensional space and xi x ) is a generic function. Such a term 
would result, for example, by imposing a constraint on the slopes of the surface with respect to a reference substrate 
plane. This is often realized in restricted solid on solid (RSOS) models for interface growth (Kim and Kosterlitz, 
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1989). These are lattice models in which the height of the surface from the reference plane at two neighboring sites 
can differ only by a given number of units, for example, ±1. The effect of such a constraint propagates at distances 
larger than the lattice spacing by discouraging configurations with large slopes (in view of their low "entropic" weight). 
Hence, at a coarse grained level, the effect of the constraint can be expressed by a term derived from Eq. ( p8| ) in a 
continuum description. The functional derivative of Eq. ([38]) is performed in detail in Appendix |e] with the result: 



Qz = — ■ 



H 



\/9 Sr 

dn z 



X( n z) 



d 2 x{n z ) 
dn z 2 



d\z-r)d inz . (39) 



The first term in Eq. ( p9| ) is the one we would obtain by assuming that d D ' s-yfgn- z, the projection of an infinitesimal 
area onto the substrate plane, does not change under the transformation r — ► r + Sr. In the Monge parametrization 
(see Eq. (jE6|)) the growth equation becomes (n z — l/^/g): 

dth = yfgQ z 

'Vh\ ( dx(n z ) \ d 2 x(n z )V_h-V 1 



V9 



y/9/ \ dn z ) dn z 2 g 



= (xi - xi)[v 2 ft + \v 2 h{vhf] + - x 'i xi'MYMY^) 2 ) + ■■■ (40) 

where the last line contains a small gradient expansion with xi = x(l))Xi = dx(x)/dx\ x= i etc. Note that even 
though the first two terms in Eq. (MG) are the same as those in Eq. (Efl) for the surface tension, the latter equation is 



rotationally invariant while the former is not, due to the presence of the other terms in Eq. (|40|). Eq. (40) also implies 
that the linear relaxation of the interface (via the term V 2 /i) is related to the behavior of x( n z) near n z = 1. From 
our previous discussion, one expects that the effect of the constraint in RSOS models is negligible on flat regions, 
i.e., Xi — 0- If is a l so reasonable to assume that x( n z) is an increasing function of n z which attains its maximum at 
n z = 1. This would indeed favour flat regions over inclined ones. This would therefore lead to a positive coefficient 
Xi - Xl > in Eq. ©. 

In principle one cannot exclude situations in which Xi ~ Xi < 0- This would lead to a linear instability in the 
growth equation of the same form as the one which has been invoked by Mazor et ai, (1988) to describe the columnar 
morphology of thin films (see later). 

It is interesting to note that in D = 1 the two nonlinear terms in Eq. ( f4(i| ) have the same form (= (dh/dx) 2 d 2 h/dx 2 ) 
and the coefficient Xl ~ x'i — §Xi may assume both positive and negative values while keeping x'i ~ Xi > 0. 
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Here and in the following, if not stated otherwise, we assume that the direction z is the same as that along which 
h is measured. The case in which a different direction is singled out can be dealt with similarly. A specific example 
will be discussed in the next section. 

7. The flux of particles and geometric effects 

Apart from forces that can be derived from an Hamiltonian, effective forces resulting from purely geometrical 
effects may also appear in Eq. (|20|). The simplest such non-conservative term derives from a flux of particles, whose 
velocity is $, that reaches the surface and sticks to it. The external flux of incoming particles is the basic source 
of non-conservative noise in interface growth. For this reason, we consider here only the effect of the average flux 



J = ($) while the fluctuation term F = $ — J will be the subject of section 1II.B . The growth rate Qj produced by J 
is proportional to the flux of J through the surface, and is exactly a measure of how many particles have been added 
to the surface: 

Qj = -n- J = Jnj, (41) 

where nj is the component of h in the — J direction, J = \ J\ and the negative sign recalls that n and J have opposite 
directions. When growth occurs from a vapor, the average of at the position f of the interface is proportional to 
the normal h at r, i.e., ($) = —An, implying C? vap = A. Thus, growth in the presence of vapor has the same effect 
as the pressure term in subsection 2. On the contrary, when particles arrive at the surface in a collimated beam this 
term can be eliminated by a Galilean transformation. Indeed, in the framework of Monge parametrization one gets: 

c^lnux = VaGj = —J z + ±±_ • v/i, (42) 

with Jj_ being the projection of J on the substrate plane. The first term is absorbed into df h by the transformation 
h — > h — J z t while the second disappears, once x — ► x + Jj_t. 

Let us return to the height constraint used in some discrete lattice models of the growth of interfaces. Apart from 
producing a surface tension term, the constraint also has an effect on the flux term. Let us consider a vertical flux 
J || z for simplicity (Jj^ = 0). At a microscopic scale, deposition can only occur in local minima of the discretized 
surface. In a coarse grained picture, minima are rare on steep portions of the interface and deposition would be less 
probable there. In other words, the constraint should reduce the flux d(f> = njdP s^fg through the infinitesimal surface 
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element dPs^/g, if nj = n z is small. This is easily accounted for by multiplying Eq. (^l|) by an increasing function of 
n z , T(n z ) 

Qj = Jn z T(n,). (43) 

In the small gradient expansion in h, this term produces the nonlinear term of the KPZ equation for any choice of 
the function T(a;), provided it is increasing. Indeed, if J || z, one has 

d t h\ flux = J^gn z T(n z ) = JT(1 - (Whf/2 + . . .) 

which, in a Taylor expansion, contains the non-linear term of the KPZ equation with a negative A. This result was 
also found by other means (Krug, 1989; Tang et at, 1992). The effect of a constraint on height differences thus 
produces both the surface tension, as seen in the previous section, and the nonlinear term of the KPZ equation. 

Of course the applicability of Eq. (^) extends to any situation where the probability for an incoming particle to 
stick on the surface depends on the local inclination of the surface. This can also account for the relation between the 
ballistic deposition model and the KPZ equation. In this model (Meakin et at, 1986), particles come from a beam 
downwards in a straight line and attach at the first site they reach in their trajectories which has a nearest neighbor 
surface site. The result of this mechanism is a non-compact cluster with a fixed density. Even though overhangs 
are present in the surface, the description of the process in terms of a single valued function h[x, t) is possible, at 
a mesoscopic level, by considering h(x,t) as the z coordinate of the highest occupied site for each x. An incoming 
particle may either stick at the surface or penetrate into the voids of the structure. In the latter case, the deposition 
process will not result in an increase of h(x, t). Particles that arrive on a flat portion of the interface have a higher 
probability of penetrating into the structure than those arriving on steep ones. This situation would also be modeled 



by Eq. (43), with a function T(n z ) that, in contrast to the previous case, should now be a decreasing function of 
n z ; in the small gradient expansion, this would finally result in a non-linear term of the KPZ type with a positive 
coefficient A. 

The flux term produces another interesting effect in situations where the atoms cannot be approximated as point 
particles. Mazor et ai, (1988) have shown that the finite size of particles plays an important role in MBE experiments 
in which thin films are grown at intermediate temperatures. The basic observation is that, if the atoms have a radius 
£, deposition does not actually occur on the surface but at a distance £ from it in the normal direction. The growth 
rate is then proportional to the flux of the beam through a surface r'(s)= r(s) +£n that is displaced by an amount £, 
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in the normal direction, from the actual surface. The growth rate on a surface element da = d D s^/g is proportional to 
the flux d(f> = —n' ■ Jyfg 7 d D s of J through the surface r ' . Here the primes refer to the displaced surface. Since n' = n, 
the only effect comes through the ^[g factor. The metric tensor g'^ is obtained by observing that c^r ' = dif+ £,dih 
so that g'^ = dif ' ■ djf ' — gij — 2£6i.j + 0(£ 2 ) (see Appendix [A|, especially Eq. flA5|)). Evaluating the determinant 
of g'ij, we find that Eq. (^Tj) has to be modified to: 

= -h-jJ^ = -h-J{l-tH)+0(?). (44) 

In the equation for h(x, t) a term: 

y /gQt = -vg^ = = -v<:V 2 h+... (45) 
V9 

must be included in the growth equation, with J || z and = which is positive. Thus, the finite size of the 
incident particles gives rise to an effective anti-diffusive behavior in the growth equation. 

An interesting extension of the results of this section arises when considering a beam of particles which is not 
perpendicular to the substrate. This setting has also been considered by Meakin (1988b) and Krug and Meakin 
(1989,1991) for ballistic deposition processes. From our perspective, one generally expects that, in this situation, a 
coupling between the interface fluctuations and the flux term arises. Let us consider, for example, how an oblique flux 
modifies the effect of the non-zero radius of particles just discussed. It is straightforward to derive the small gradient 
expansion of Eq. ( fti| ) in the case J = (Jj_, J z ), where the z direction is normal to the substrate: 

y/gGt = Jz - J ± • Vh - i3,y 2 h + £ J ± ■ VhV 2 h + ... (46) 

As already mentioned, the first two terms are eliminated by an appropriate choice of the reference frame. The third 
term in Eq. (^) has just been discussed. The last one is new and it represents the coupling of interface fluctuations 
with the transverse component of the oblique flux. It has been argued by Marsili et al., (1996), applying dynamical 
renormalization group techniques (see, e.g., Ma and Mazenko, 1975; Forster et al., 1977; Frei and Tauber 1994) that 
this term, in competition with surface diffusion, is responsible for a new scaling behavior of interface fluctuations at 
intermediate scales. This approach predicts a roughness exponent a — 1/3 which is in very good agreement with 
experiments of MBE, where a value a — 0.30 — 0.33 was measured (Salvarezza et at, 1992; Herrasti et at, 1992). 
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One can, in general, imagine other effects arising from an inclined flux. For example, if J_j_ / a new term would 
also arise from the effect of a constraint on height differences on the flux term described in Eq. (f43|). Indeed, if 
n ■ J = n z + n_i_ ■ Jj_> one would expect a term ( • V/i)(V/i) 2 in the gradient expansion. Common wisdom ( i.e., 
power counting) however, suggests that this term will not modify the leading scaling behavior of interface fluctuations 
in the presence of the KPZ non-linearity (V/i) 2 . 

8. Surface diffusion 

In cases in which the binding energy of particles on the surface is large compared to thermal energy fluctuations, 
the motion of the particles is constrained to be along the surface. This is actually the case in many experimental 
realizations of MBE. A force acting in a direction normal to the surface cannot displace the particles. Its effect is to 
change the local chemical potential fi. Differences in chemical potential on the surface, in turn, produce a current 
proportional to the gradient of p, on the surface (Villain 1991). The evolution of the surface is governed by the 
continuity equation relating the particle density p and this current: dtp oc n ■ dtf yields the left-hand side of Eq. j2C| ) 
while the divergence of the current gives the right-hand side. 

The mathematical translation of this argument is as follows. The constraint that the motion of particles occurs on 
the surface implies that the volume V it encloses cannot change. Since 



d t V = I d^s^—jdtris, t) = I d u s^h ■ fiG + F{x, t) , (47) 



the condition dtV = for the deterministic part of Eq. (20) implies: 



d t V = J d D s^g = 0. (48) 

A sufficient condition for this to hold is: 

Q = - AT= -divj. (49) 

In other words volume conservation requires Q to be the (covariant) divergence of a (contravariant) current J % ( see 
Eq. |A3| ). Here J_ is the surface current of particles and T is proportional to the chemical potential. This relaxation 
mechanism is known as surface diffusion. If this force, JF, derives from a potential, then: 

n-T= - — . (50) 

y/g or 

22 



Having derived, in the previous section, the forces coming from the simpler contributions to the potential energy of 
an interface, it is straightforward to find the appropriate Q° term in the equation for d t f under surface diffusion by 
applying the Beltrami Laplace operator to the terms derived previously. In the Monge form, the term appearing in 
Eq. @ has the for m: 

VgG c = -Vg^ = VaA 6 -^ (5i) 

where the second equality holds for growth mechanisms that can be derived from a potential (using the result of 
Appendix |b|). When surface diffusion occurs to minimize the surface area, the corresponding term in Eq. (20) reads: 



G c s = -VsA [h-Ar )=-VsAH. (52) 



This follows from the two previous equations and Eq. (|25|). This term has been widely used in numerical and analytical 
studies (Mazor et al, 1988; Golubovic and Karunasiri, 1991; Siegert and Plischke, 1992; Sun and Plischke, 1993; Sun 
and Plischke, 1994). Its expression, in terms of h(x,t), is 

V5^ c = -MsV5AV-=^~^ s (V 2 ) 2 ^ + ... (53) 

A pressure term does not produce any surface diffusion, since surface diffusion entails the conservation of the 
volume. External potentials, that depend on r, can produce surface diffusion. In the simplest case of the gravitational 
potential (along the z direction) we have: 

Q% = v g z ■ n H (54) 

where v g = pa g is the mass density times the gravitational acceleration. The explicit n ■ z factor indicates that this 
term breaks rotational invariance. Q g has the same form as in Eqs. ([11]) and (^5|), but has the opposite sign. In 
actual experimental situations however, the gravitational energy is negligible with respect to the binding energy. It 
can be estimated, for example, that in MBE growth of thin films v g is of the order of 10 -14 of the corresponding 
coefficient in the term. However, terms of the form: 

V9Gg = VgV9Ah = u g V_ ■ =- ~ u g V 2 h +... (55) 

with a positive v g have often been used in recent publications (see, e.g., Golubovic and Karunasiri, 1991). This is 
principally justified by renormalization group considerations. Even if this term is not present at a microscopic level, 
it is generated in the iteration of the renormalization group equations from the non-linear terms. 
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Surface diffusion can also be induced by a curvature dependent Hamiltonian or by orientation dependent potentials 
in the same manner. The explicit expression is simply given by that for the non-conserved case with the additional 
Bcltrami-Laplacc operator. 

B. Stochastic Evolution 

As already mentioned, the principal source of randomness in interface growth comes from a flux of particles that 
deposit on the surface. Another kind of noise is produced by thermal fluctuations of the surface. The main difference 
is that the latter conserves the total volume enclosed by the surface. As discussed previously, the random force 
F = nrj is in the normal direction and we can take (77) = 0. In some representations, e.g., if one studies the growth of a 
spherical substrate in polar coordinates, this choice cannot be made. The properties of 77 will now be discussed for the 
two different cases of non-conservative and conservative noise. Our main concern is the correlator (j](s,t)r](s',t')}. 
In cases when the statistics of the noise is Gaussian, as is almost always assumed, this correlator specifies the entire 
distribution of rj{s,t). (Non-Gaussian statistics for the noise, especially with long tails, has been shown by Zhang 
(1990) and Krug (1991) to affect the scaling properties of the interface). 

1. Non-conservative noise 

For a flux of particles arriving at the surface with velocity the noise term is given by 77 = n ■ F where F = $ — J, 
with J = (<E>). It describes the fluctuation in the number of particles arriving at the surface. The correlator of 77 is 
then given by: 

(r,(s,tMs',t')) = n a {s,t)n (j {s\t'){F a {lt)F^{s',t')) 

= njs, t)na(S l!)T° ^ a ~J ' Sit - if). (56) 

Here we choose F to be delta correlated both in space and time. The more general choice: 

(F a (s,t)F^(s',t')) = T afj C{s-s\t-t'), 

where C(s,t) is a R-invariant function, involves no further complication, so we will use the explicit expression (|5^). 
The coefficients T a @ are symmetric in the indices and specify the geometric properties of the noise. We may distinguish 
two extreme cases: 

• deposition processes that occur from the condensation of an isotropic vapor and 
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• deposition that occurs from a collimated beam of particles. 

The difference between these two possibilities fully appears in the Monge parametrization. In order to discuss this, it 
is useful to introduce the random field: 

f){x,t) = y/gr](x,t), (57) 

which is the quantity that appears explicitly in the equation for d t h. The correlation properties of f){x,t) are derived 
directly from those of rj(x, t) discussed above. 
^From equation (|5^) one easily sees that: 

(r](x, t)f](x ,t )) = — o(x — x )o(t — i J. (58) 

If growth occurs from the condensation of an isotropic vapor, we expect that F is a random vector with uncorrelated 
components andQ T Q/3 = T5 a/3 . In equation (|58|) we find: 



(fj{x, t)fj(x!, t')) = T^8(x - x!)6(t - t'). (59) 

For growth occurring from a directed flux, one may assume that all components of F are independent random 
variables so that T a ^ = T a S a ^ (no summation on a is assumed here). If rotational invariance is expected for 
rotations in the substrate plane wc have: Ti — Ty and 

{nix, t)fj(x', t')) = Tz+T ^- K)2 5{x - x') S(t - t'). (60) 
V9 

Note that in-plane correlations are enhanced in regions where h has steep derivatives. For a collimated beam per- 
pendicular to the interface we have T\\ <C V z and 

(fj(x, t)fj(x!, t')) S ^pS(x - d) $(t - 0> ( 61 ) 

while if T|| = T z we recover Eq. ([59]). If the randomness only affects the intensity of the beam, T a ^ — TJ a J^ and 
equation (^8|) becomes: 

(fj(x, t)fj(x!, t')) = r (Jz ~-^'- ±) 5{x - rf) S(t - t') (62) 

V9 



1 Here we use the Kronecker symbol with both upper indices which coincides with the metric tensor g al3 in the D+l dimensional 
space. 
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where Jj_ is the component of J in the substrate plane r — x. In the case of vertical rain, J z 3> | </j_ | and we recover 
the previous result Eq. 

The physical meaning of the prefactors of the delta functions is evident if we introduce rj {x,t) such that 
(rj (x, t)rj (x^ , t')) = T6(x — x^)S(t — t'). In condensation from a vapor, Eq. (]59|), we find that ?j{x, t) — g 1 ^ A rj (x_, t) in 
Eq. ( pl| ) (Since g contains stochastic variables, this is also the correct mathematical interpretation of the correlation 
in Eq. (p9|)). The noise is enhanced in regions where h has steep derivatives since the exposed surface area in the 
substrate element dPx is larger by a factor of ^fg. The opposite case is that of growth from a perpendicular beam 
J = J z 'z in which case fj(x,t) — <?~ 1//4 f?o(a;, £)■ This is because the flux of J through the surface is proportional to 
h ■ z — 1/^fg and regions with high slopes receive less particles than those that are natter. 

2. Conservative noise 

Another source of noise comes from thermal fluctuations and from internal degrees of freedom of the interface. In 
this case the noise is called conservative because it causes no increase of the volume enclosed by the interface. Using 
Eq. ([47|), this requirement is translated into the condition: 

dtV\ nmse = J d D s ^ v = (63) 

where again we have taken F — nrj. This poses a condition on 77. A general way to let the noise contribution in Eq. 
( |63| ) vanish is to take: 

7} = diy C, 

where div is the covariant divergence acting on the vector £ that is a delta correlated noise both in space and time: 

<Ci(s > *X J 'Cs , ,0) = ra/fcf^ - 1'). 



Here reparametrization invariance has been satisfied and the delta function allows the use of g' = g(s') instead of 
g(s). The correlations of rj readily follows: 

(v(s,t) V {s\t')) = J=d,:{v3d; U', *'))]} 

= r-L=a, ,g<iu)'Ms - i)5(t - 1')] 

V99 

V9 



-r 

-TA 6 -^p-5(t-t') (64) 
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where the primed quantities refer to s' and the presence of the delta function has been used repeatedly to change 
from primed to unprimed quantities (note that the operator A does not act on g'). This is the natural generalization 
in reparametrization invariant form of the correlator often used in dealing with conserved noise that contains the 
Laplacian operator acting on the delta function. The expression of the correlator in the Monge parametrization is 
readily derived from the above expression. 

3. Approach to equilibrium 

In previous sections we have frequently dealt with terms in the deterministic part of the growth equation of the 
form: 

M h(x,t) = -V 9 T IK( ^ + fj(x,t), (65) 

where T = 1 for non-conservative dynamics, and T = — A for conservative dynamics. If the noise is Gaussian with 
correlation (see Eq. (|59|)): 

(r)(x, t)rj<£, 0) = 2T^gTS{x - x')5(t - t'), (66) 



one can write the associated Fokker-Planck equation in the form: 



T- 6 m 



Sh(x, t) Sh(x, t) 



P[h(x),t}. (67) 



Eq. (|67|) holds if a suitable regularization has been chosen such that Sh * x ^ commutes with y/gT. This is fulfilled, 
e. g., in the dimensional regularization scheme commonly used in field theoretical treatments (see, e.g., Zinn- Justin, 
1993). 

In Eq. (|67]) P[h(x),t] is the probability functional which yields the probability of the interface configuration h(x) 
at time t. It is easy to see that the stationary distribution, obtained by setting the right-hand side to zero, is given 

by 

P[h(x), t ->■ oo] ex exp j - ^ ( ~ )] j . (68) 

The equality of the coefficient in front of the functional derivative and of that in front of the noise correlation is often 
referred to as the fluctuation dissipation theorem (see, e.g, Deker and Haake, 1975). When this theorem holds the 
stationary distribution is given by Eq. (|68|). It is interesting to note that the R-invariant form of the functional 
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derivative, with the R-invariant form of the correlator deriving from an isotropic growth mechanism (T = 1) yields 
a dynamics which leads to the equilibrium distribution Eq. (|6^). This was first noted by Bausch et a/., (1981). One 
may apply these considerations to the R-invariant form of the KPZ equation describing growth from condensation of 
a vapor as discussed by Maritan et at, (1992). However, as already pointed out, TIkpz of Eq. ( |3l| ) is unbounded as 
h — > oo and therefore Eq. ([68]), not being normalizable, is meaningless. 

Note also that if the flux term breaks rotational invariance, as in the case of a collimated vertical beam (see Eq. 
(|62|)), the fluctuation dissipation relation does not hold even if the deterministic part derives from a potential. 



Furthermore we note that, for T = — A, the deterministic part has the same form discussed in section III. A. 8 and 



the correlation of the noise is exactly that derived above in Eq. (64). Therefore we can conclude that a deterministic 
conservative dynamics with a conservative noise leads to the stationary state described by Eq. ( |68| ) (provided P[h,t) 
is normalizable). 

IV. DISCUSSION 

We have seen that there are four different mechanisms producing a Laplacian term i/V 2 /i in the small gradient 
expansion of Eq.(pl|). Three of them, the surface tension, an orientation dependent potential (or a constraint on n z ) 
and surface diffusion induced by gravity, lead to a positive v coefficient that drives the evolution towards flatter and 
flatter surfaces. Note that one can distinguish between these effects only through higher order terms in the gradient 
expansion, which may however be irrelevant in the renormalization group sense. The fourth mechanism, related to 
the finite size of the aggregating particles, gives a negative contribution to the coefficient in front of the Laplacian, 
so it would produce an instability if acting alone. This is evident since Eq. ( pl5| ) can formally be derived as a surface 
diffusion induced by a negative, gravitational field [y g — — z/j in Eq. (j55[)). This also implies that this term strictly 
conserves the volume enclosed by the surface (while surface tension, Eq. (^6|), does not) even though it was not 
derived from conservation considerations. 

Secondly we note that the non-linear term of the KPZ equation (|l7|), can be derived in one of three ways: from 
a pressure term in a potential, that gives a positive A for a growing surface, from growth due to condensation of a 
vapor or from an inclination dependent factor in the flux term. The latter may result from the effect of a constraint 
on n z , yielding a negative A, in agreement with known results on restricted solid on solid models (Meakin, 1993). It 
has also been argued that such a term is expected in ballistic deposition but with the opposite effect, i.e., inhibiting 
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growth on flat portions of the surface. A positive A is expected in this case. The change of sign in this coefficient does 
not change the character of the process dramatically as it does for v. The value of A however is a directly measurable 
quantity (Krug, 1989) since it is related to the inclination dependence of the average velocity of growth. In this way 
Krug (1989) was able to predict the presence of the non-linear term of the KPZ equation in various models. Note 
that this is a criterion based on the global behavior, while our analysis is based on the local properties of the growth 
process. Another derivation of the KPZ equation for restricted solid on solid models is based on its relation with the 
directed polymer problem (Kardar and Zhang, 1987, Fisher and Huse, 1991)in a random environment as shown, e.g., 
by Tang et al, (1992). 

Equation ( |36| ) provides a further physical derivation of a term (V 2 ) 2 /i in the equation for a growing interface. This 
term has usually been associated with surface diffusion. This extends the validity of the results derived in the presence 
of this term to situations where the restoring force is derived from a potential corresponding to the surface energy 
being proportional to H 2 . The expansion of the potential energy of the interface in powers of H will, in general, 
contain also a constant term H , proportional to the surface area. Straightforward analysis shows that in the presence 
of this term all higher powers are irrelevant from a renormalization group point of view. There are situations, as in 
the dynamics of fluid membranes, where this term is known to be absent. The leading terms would then result from 
the terms linear and quadratic in H. An interesting point for future investigation is an analysis of the effects of the 
term resulting from a potential linear in H , Eq. (|34|). 

Sun et al., (1989) have studied an equation with a linear term proportional to (V 2 ) 2 /i and a non-linear term 
proportional to V 2 (V/i) 2 . The same equation was studied by Lai and Das Sarma (1991) and Kim and Das Sarma 
(1994) (see also Das Sarma et al., 1996). This is in a loose sense the conserved version of the KPZ equation. We 
note here that while the first term could be derived from surface diffusion or from a potential proportional to H 2 , 
the non-linear one is not related to one of the simple mechanisms discussed here. In particular, a derivation of the 
conserved KPZ equation cannot follow the same lines described above. In one case, the KPZ equation comes out as a 
result of a process that does not conserve the volume. In the second, it derives from the effect of a constraint on height 
gradients and one needs to motivate the rather odd choice T(n z ) = — AT(n z ). There are no R-invariant potentials 
that would lead to such a term in the gradient expansion, either for conservative or non-conservative dynamics. It is 
still possible however that such a term is generated dynamically in a renormalization group procedure. 

Note also that in the noise term, there are prefactors in front of the random field. In a small gradient expansion 
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these yield a nonlinear term of the KPZ type with a random valued A. It is also important to note that the occurrence 
of multiplicative noise could significantly change the scaling properties of growing interfaces. 

We also discussed the generalization of the fluctuation dissipation theorem to the R-invariant form of the growth 
equation. We recovered the observation of Bausch et at, (1981) that R-invariance and rotational invariance are 
enough to establish fluctuation dissipation relations for non-conserved dynamics. We also found that this observation 
generalizes to conserved dynamics (both the noise and the deterministic part being conservative). 

Recently, Keblinski, Maritan, Toigo, Koplik and Banavar (1994) and Keblinski et at, (1996) have introduced a 
simple continuum model that allows for overhangs and an arbitrary topology of the growing interface. The model 
captures surface diffusion in a natural manner and, with an appropriate aggregation mechanism, it produces growth 
normal to the interface. The model equation consists of two parts, the first is conserved order parameter dynamics 
that allows for the definition of a topologically unrestricted interface and builds the correct physics of surface diffusion, 
whereas the second term provides the growth and roughening at the interface. 

In the simplest version, their equations are: 

9mt) -TV>J^ + I, (69) 



dt Sf(r,t) 



where 



F = 



-If + ll+aiVff 



dv, (70) 



and I = C 1 \Vf\+D 1 ^/\Vf\ri(r,t), (71) 

where r\ is a Gaussian noise uncorrelated in time and space with a width equal to 1 and mean value 0. 

Equation ( |69| ) without the / term has a simple interpretation - it is merely the deterministic part of the standard 
model - B dynamics (Hohenberg and Halperin, 1977) that conserves the order parameter. 

The choice of the sign of the coefficient of / 2 in the expansion for the free energy (^) corresponds to a temperature 
lower than T c , so that the two values ±1 of the order parameter f{r,t) minimizing the free energy describe the two 
equilibrium phases of the system. 

An interface can naturally be defined as the crossover region between the / = — 1 and / = +1 regions - and, 
operationally, a point fi is defined to be on the interface when f(fi,t) = 0. 

a > is the surface diffusion coefficient and sets the intrinsic length in the system. Indeed the width of the interface 
is proportional to y^a. Also the effective strength of the surface tension turns out to be proportional to ^fa. 
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The I term allows for the growth and fluctuation of the interface. The V/ factor ensures that the growth and 
fluctuations are operative only in the vicinity of the interface - away from the interface V/ is effectively zero. The 
positive coefficients C\ and D\ are the magnitudes of the growth and noise respectively. The V/ factor produces 
growth normal to the interface. 

Numerical results (Keblinski, Maritan, Toigo, Koplik and Banavar, 1994; Keblinski et al. 1994, 1995, 1996) show 
that this model is in the same universality class as the KPZ equation. The / growth mechanism not only gives 
growth normal to the interface but also the rate of the growth per unit length is constant along the interface, equal 
to C\ J!!™"* \Vf\ds — 2C\ where the integral is performed across the interface and s m i n and s max are defined by 
f{s m in) = — 1 and ,f{s max ) = 1 (note that the integral is independent of the shape of the V/ profile as long as / 
increases monotonically from / = — 1 to / = +1). This feature leads to KPZ-like behavior. 

The model presented above has a local conservation law to avoid the formation of islands with / = — 1 in the 
region predominantly with / = +1, and vice- versa. However, the conserved and non-conserved model should exhibit 
interfaces with the same behavior. It can be shown (Keblinski et at, 1996) that the non-conserved version of this 
model is equivalent to Eq. ( p0| ) with a noise term of the form ( |57j ) with variance given by ( |59| ) when the interface is 
sharp, i.e., a is small. Thus, this model may be interpreted as a continuum version of the Eden growth model, with 
redistribution of the aggregated particles via surface diffusion. The Eden (1958) model is known to be in the KPZ 
universality class - the surface diffusion in the limit of large length scales does not change the geometrical properties 
of the interface, but introduces a short range smoothing mechanism. 

The growth of real surfaces is often influenced by non-local effects such as screening or shadowing. When the 
aggregating particles follow linear trajectories one can expect that, if the roughness is large enough, some parts 
of the interface are shadowed and, therefore, do not grow (Tang and Liang, 1993). In order to accommodate this 
phenomenon, the previous model can be extended (Keblinski et al, 1995, 1996) to incorporate the dynamics of the 
depositing vapor and non-local effects. 

The extended model involves two fields / and 4> and is governed by the equations: 



+CVWJNv(r,t) (72) 
= V[DV<f>(?t) - Mr, t)} - B(V/) 2 ^(f, t) (73) 



dt 
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with F again given by Eq. (|70[). While the first part of Eq. ( |72| ) is identical to that described earlier, the growth 
mechanism is different. Now the growth of the / field occurs at the expense of the <fi field. The <f> field represents the 



local density of the incoming particles towards the interface, and Eq. (73) describes the dynamics of the depositing 



vapor. The first part of Eq. ( |73| ) is simply the diffusion equation in the presence of an external force A. In order 
to analyze the growth arising from ballistic trajectories, D has to be chosen much smaller than A, so the A<f> flux is 
the primary mechanism for 4> field transport!!. The aggregation relies on the conversion of the (f> field into the / field 
as described by the coupling term B in Eqs. (^) and (^3|). The V/ factor in the B term makes the aggregation 
operative only within an interfacial " skin" region of the aggregate with its width proportional to y/a. The <f> factor in 
the aggregation term ensures that the growth occurs only if cf> > 0. The B term acts as a sink for the diffusive field 4>, 
and its magnitude is chosen to be sufficiently large to convert all of the (j) into / within the interfacial region effectively 
leading to ~ below the interface (/ ~ +1). Shadowing effects are naturally incorporated in the equations. When 
the 4> field trajectory intercepts the skin, the <f> field is converted into / and any subsequent interception occurs with 
= and, therefore, does not lead to the growth of the shadowed part of the interface. Note that non-local effects 
are incorporated in a local way in Eqs. (|7^) and (|7^). One does not need to monitor the geometry of the interface to 
incorporate shadowing - it is implemented dynamically by the field. The conversion of the field into the / field at 
the vicinity of the interface does not depend crucially on the particular functional form of the coupling chosen. The 
rate of growth is effectively equal to the intensity of the incoming field flux. In addition to the growth term, there 
is a fluctuation term C in Eq. ([72]). The Gaussian T)(f,t) factor is the same as was introduced previously. In this 
manner, the fluctuations in the strength of the incoming flux are incorporated, since the aggregation rate is equal 
to the intensity of the incoming flux. 

Strikingly, the model presented in Eqs. ( |72|) and ( f73[) can be straightforwardly modified (Keblinski et at, 1994, 
1995, 1996) to model diffusion-limited-aggregation type of phenomena (Witten and Sander, 1981; Niemeyer et al. 
1984; see also Pietronero and Tosatti, 1985). One just needs to substitute the ballistic flux with a diffusive flux that 
is responsible for the transport of the aggregating field. 

The equations are modified to: 

^M-V 2 ^- + J(rt) (74) 



2 Here D is not to be confused with the spatial dimension. 
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d<f>(r,t) _ r ^ 2 



at 

with 



= DV"<f>(?,t)-J(r,t), (75) 



J(r,t) = -Vf ■ DV(f>-f](r,t), (76) 

where rj{f, t) is a Gaussian noise with non-zero average value V > and width W. 

The first part of Eq. ( f74| ) and the free energy F is the same as described previously. The interaction term J > 0, 
leads to the growth of / and decay of cf> such that / + (f> is a conserved quantity and changes only due to the sources of 
the <p field at the boundary. The R-invariant form of the growth equation provides an interesting alternative approach 
to avoid the no-overhang approximation, even though it cannot describe non-local effects like shadowing in ballistic 
aggregation and screening in the diffusion limited aggregation models. 

The inadequacy of the no-overhang approximation is just one of the reasons why one may want to rely on the 
R-invariant formulation of growth problems. For example, the R-invariant form of the equation displays the full 
invariance properties with respect to space translations and rotations and the conservation laws which the process 
satisfies. These are lost, for any but infinitesimal transformations, once one restricts attention to the lowest order 
terms in the gradient expansion. These invariances play a crucial role in the implementation of renormalization group 
approaches around the lower critical dimension (see, e.g., Bausch et a/., (1981), because they provide conditions for 
the renormalizability of the theory. 

A second situation in which a full R-invariant form of the growth equation would be preferable to the lowest 
order gradient expansion arises when the scaling behavior is determined by a strong coupling fixed point. In order 
to appreciate this situation it is preferable to sketch briefly the standard approach of the perturbative dynamical 
renormalization group (Ma and Mazenko, 1975). Let us consider the Langevin equation 

d t h(x, t) = C h(x, t) + XAf[h(x, t)] + r)(x, t) 

where £ is a linear differential operator in the x variable and J\T is a non-linear combination of h and gradients. As 
an example, in the KPZ equation ( |l7| ) one has £ = V 2 and Af[h] = (V/i) 2 . The noise is Gaussian with zero mean and 
(r](x,t)r](x! ,t)) = 2T5 D (x — x')S(t — t'). If A = 0, the equation can easily be solved in Fourier space. It is therefore 
easy to find the exponents ao and zq which characterize the scaling in the linear theory: Under a rescaling of length 
by a factor of £, the time will scale by a factor £ z ° and h will acquire a factor £ a °. One can then analyze the effect of 
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a small non-linearity (A <C 1) on the dynamics. If the above discussed change of scale in the linear theory affects the 
non-linearity with a factor l v ( i.e., if Af[h] — > £ v Af[h] as x — > %, < — > ^ z °t and ft- — * £ Q °/i) and y < 0, one can conclude 
that the nonlinearity is irrelevant. On the other hand, if y > one concludes that the effect of the non-linearity 
will increase as the scale increases and will eventually dominate the large scale behavior of interface fluctuations. It 
usually happens that y is a decreasing function of D, the dimensionality of the substrate, and there is a dimension 
D c below which the non-linearity is relevant. In this case, the full program of the dynamic renormalization group 
becomes necessary. Since the method can in principle treat only small non-linearities, being based on a perturbation 
expansion, it usually provides estimates of the exponents close to the dimension D C1 where one can assume that since 
y is small, the fixed point is accessible within the expansion. If this procedure works, i.e., if one finds a stable fixed 
point whose distance from the fixed point of the linear theory is small when D c — D is small, one can also conclude 
that all higher order terms, which have been neglected in the gradient expansion, are irrelevant (if ao < 1). This 
whole program, which is a very powerful tool to estimate critical exponents, fails if one finds no stable fixed point or 
if the non-linearity turns out to be relevant only if A is bigger than a critical value A c . The nature of the phase for 
A > A c turns out to be outside the range of perturbative methods and in this situation, one has no reason to neglect 
the higher order terms in the gradient expansion. 

This situation is realized in the most studied model of growth, namely the KPZ equation, for which above the 
substrate dimension D c = 2 (which is also the physical dimension in which one is interested) the non-linearity 
becomes relevant only if A > A c . 

There is a general consensus that different models for non-equilibrium dynamics, such as restricted solid on solid 
models (Meakin, 1993), the Eden model (Eden, 1958) and directed polymers in random media (Halpin-Healy and 
Zhang, 1995), fall in the same universality class. Such an expectation is mainly based on the observation that the 
small gradient expansion of these models contains the terms in the KPZ equation ([T^). While analytical and numerical 
results in D = 1 almost unambiguously support this expectation, there is no reason, in principle, to believe that in 
D > 2 the strong coupling regime of these models is described by the same exponents. In fact, as we have shown, 
the continuum equations for the Eden model, for RSOS models and for ballistic aggregation differ in the higher order 
terms of the gradient expansion. A yet different gradient expansion can be obtained for the free energy of directed 
polymers. Going one step further than the usual KPZ truncation, Marsili and Bray (1996) studied the equation 
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21 = v ^h + n{Vh) 2 V 2 h + A + - (V/i) 2 + . . . + r]. (77) 



This derives from the small gradient expansion in the case of RSOS growth, as described in section III. A. 6. This 
equation indeed coincides with Eq. ( |40| ) with \i = 0, xi = ~Xi > 0- Note that here both surface diffusion and the 
growth term have been expanded to the same (second) order. It turns out that a mean field, infinite dimensional 
limit of Eq. @ can be meaningful only for k > (Marsili and Bray, 1996). This analysis reveals that steps or 
bumps of a finite height develop on the surface and that dynamical scaling (Family and Vicsek, 1991) is not satisfied. 
These results, which differ substantially from those obtained in the same limit for directed polymers by Derrida and 
Spohn (1988) (see also Cook and Derrida, 1991), raise doubts on the existence of a single universality class for all 
these processes. In high dimension, the k term turns out to be necessary in order to avoid finite time singularities 
which occurs in the simpler KPZ equation (k — 0). These finite time singularities turn even worse in the case of the 
Eden model or ballistic aggregation. Indeed, the model Eq. ([77]) can be generalized to these process with k < which 
describes the suppression of surface tension on steep portions of the interface discussed by Maritan et al, (1992). 
These instabilities suggest that the Monge representation of the interface may not be adequate to describe these 
processes (Maritan et at, 1992). 

Anomalous dynamic scaling has also been reported ( Das Sarma et al., 1996) for various models proposed to describe 
MBE. More precisely, numerical results show that the local scaling properties, defined in terms of the correlation 
function, identify an exponent ai oc which happens to be different from the one which describes the behavior of the 
global surface thickness W(L, t) with the size L of the sample at saturation. This anomalous behavior has been 
related by Schroeder et al. (1993), Das Sarma et al. (1994) and Krug (1994) to the peculiar statistics of steps ( i.e., 
height differences between neighboring points). The broadness of this distribution (which diverges for infinite times 
and infinite L) again raises doubts on the validity of the small gradient expansion. 

In all these situations, we believe the concept of reparametrization invariance may prove to be an invaluable starting 
point for elucidating the correct physics. 
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APPENDIX A: DIFFERENTIAL GEOMETRY 



An orthonormal basis is assumed in D + 1 dimensional space and Greek letters are used for the vector components. 
Latin letters used as an index refer to the components of vectors in the D dimensional parametrization space. s l are 
general curvilinear coordinates that label points on the D dimensional surface. The notation di = d/ds l is used for 
covariant derivatives i. Summation over repeated indices is always assumed. Lastly, for the scalar product in both 
spaces a dot is used while x denotes the vector product. 

The distance between infinitesimally close points on the surface is given by the first fundamental quadratic form: 

\dr\ 2 = dirds* ■ d 3 rds j = g ij ds i ds j . (Al) 

This defines the metric tensor: g^ — di? ■ dj?. g = det{gij} denotes its determinant while g^ is the inverse: @ 
9ikg k ^ = The metric tensor and its inverse are also used in the usual way to lower and raise indices, i.e, V{ — gijV^ 
or v % — g^Vj. We also use the notation d l — g^dj. 

The only restriction on the choice of the parametrization is that g ^ 0, i.e., that g+j is invertible, and this implies 
also that di? ^ 0. The vectors di? lie in the tangent hyperplane so that the normal versor is given by n — g~ 1 / 2 d\? x 
82? x ... x do?, where g -1 / 2 ensures normalization (v±x ... x vn = £ ao a 1 ...a D vi ai ■■■VDa D , where £ ao a 1 ...u D is 
the completely antisymmetric Levi-Civita tensor and is equal to (— l) p with P being the order of permutation of 
ao«i ■ ■ ■ ot£, with respect to 1, 2, . . . , D + 1.) 

A quantity X^.. 3 '"(s) is said to be a tensor if, under the change of parametrization s'(s), it transforms as: 

Tl J (s') = ^...^T k J-(s). 



^From Eq. (Al) one sees that and its inverse g %3 are tensors. 

A quantity c/?(s) is said to be a scalar if <p'(s') = (p(s) . In particular r a (s) and n a (s), with a = 1, 2, . . . , D + 1 are 
scalar quantities, while diip and d z ip are particular cases of tensors called co-vector and vector respectively. 

The invariant surface element is given by da = d D s^fg and this implies that the invariant form of the delta function 
in parameter space is: 



3 We have avoided introducing covariant derivatives in order to maintain the exposition as simple as possible. Thus di and d' 
introduced here behave like tensors only when applied to scalar quantities. 

4 The symbol <5 is used here for the Kronecker delta. It will also be used for the Dirac delta function and for functional 
differentiation. 
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*oU-l') = ^M^ (A2) 

where 5(s) is the usual delta function in D dimensional space. Thus J d<rf(s)So{s — s') = f(s') 

For differential calculus, invariant forms of the gradient, divergence and curl are obtained requiring the transforma- 
tion properties of tensors to apply. The gradient of a scalar S is simply given by diS while the divergence of a vector 
isi: 

divu = ^S i (Vff« i ). (A3) 
V9 

Taking the divergence of the contravariant gradient yields the reparametrization invariant generalization of the 
Laplacian operator in curved space: 

A = -jfiy/gP = ^ ( \ Tin' ' 'K ) (A4) 
which is known as the Beltrami-Laplace operator.! 



The curvature K of the surface along a curve s(£) is given by h ■ 9 ( I is the arc length). Since n ■ dir = 0, n 



may be written in terms of the second fundamental quadratic form: 

, ds l ds j . . 

K = b ^^l' (A5) 

where bij — h ■ didjf= —dih ■ djf% This defines the principal curvatures (directions) as the eigenvalues Xi (vectors) 

of b\ . These are invariant under reparametrization. The mean curvature H is the sum of these and thus equals i the 

trace of b\ : 

D 

H = bl = J2 A * = -di™ ■ 9^- ( A6 ) 

i—l 

Another useful definition of H comes from observing that, since h ■ d i f = 0, di y-s/gn ■ <9V) = ^/g{dih) ■ d % r + h ■ 
di (y/gd l f) = 0. This implies: 

H = -d i h-d i r = n- Ar. (A7) 



5 Indeed, if ip and if are a scalar and a contravariant field respectively, Eq. 
— J dP s^/qcp div v. 

can easily be deduced by requiring J dadiipdjipg*-* = — J da<p A </> f° r an y scalar field ip. 

7 The first form shows that b is a symmetric matrix, whereas the second shows that it is a tensor: hence the R-invariance of 
its eigenvalues. 

8 Actually the mean curvature should contain a factor 1/D that we disregard for convenience. 



A3) follows on requiring J dP s^/gv 'di<p = 
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The Gaussian curvature is defined as K = det{6^} = J\ i A,-. Furthermore it is easy to see that 

dS = -bijdir. (A8) 

Indeed dih-Ln and diifi ■ djf) = implies dih ■ djf = — n ■ didjf = — £>y. 

1. The Monge Form 

A particular choice of parametrization is the Monge form: 

f>=(x,h(x)), (A9) 

where x is a vector in the D dimensional substrate plane and h(x) is the height of the surface in the direction z 
perpendicular to this plane. Use of this parametrization implies that no overhangs are present in the surface since 
otherwise h(x) would not be single valued. In this parametrization the metric tensor has the form: 

9ij — ^ij + dihdjh and g 1 ^ — rJy — —dihdjh, (A10) 

where: 

9 = 1 + {Yhf and ; h = —{-Vh, 1) (All) 

and 

I'j //"',/... <V,^4 (A12) 

Finally the mean curvature is given by 

H = bU=sZM. (A13) 
y/9 

The equation for h(x,t) is obtained from ( po|) considering the various components of r. On defining r(s,t) = 
(x(s, t), /i(s, i)), we get: 

d t h(s, t) = n z G, 
d t x i {s,t)=n i Q, 

where ri 1 and n z are the components of the normal in the directions x l and z respectively. These derivatives are 
evaluated at constant s while we are interested in the derivative of h at constant x, 
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dh 

8 t h{s,t) = d t h(x,t) + —d t x l (s,t). (A14) 

where h(x,t) = h(s(x,t),t). ^From the above equations and h = (— V/i, l)/y/g, one readily finds the deterministic 
part of Eq. (§l]). 

APPENDIX B: EQUATIONS DERIVED FROM A POTENTIAL IN THE MONGE REPRESENTATION 

The property that the functional derivative is orthogonal to the vector dif translates into: 

^—d h + — - 

Sh 1 Sx l 

This allows one to eliminate the functional derivative w.r.t. x l in 

, SH Z SH t SH ^SH 
-^g = n. — =n— + n— i =V9 W 

and to find Eq. (||). 

APPENDIX C: DERIVATION OF THE GROWTH TERM DUE TO SURFACE ENERGY 

In the functional derivative of g w.r.t. r, in Eq. (p^), we use the fact that g is a determinant and the property: 

S In dct M = Str In M = tiM~ 1 SM (CI) 

that holds for variations of a matrix M. This allows one to write: 

S J d D s^ = i J d"s x 7}g"6g„ =- J d D sd l (^d J r) ■*?=-[ d D s^gAr- Sr 

for a variation Sr of r , which readily yields: 
1 S 



. / d D s'JgJj) = -±=d i (^gJT)g^m) ?(s) — A f(s), (C2) 



where Eq. (A4) has been used. 



APPENDIX D: CURVATURE DEPENDENT POTENTIAL 



For a variation Sr in r, h changes to h + Sn. Since S(h ■ h) — 0, the variation Sh is normal to h. Consider Eq. ( A7) 
for H . In 

SH = Sh ■ Ar + h ■ S(Ar) 
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the first term on the right-hand side vanishes since Ar is parallel to n. The variation of Ar, by simple arithmetic, is: 

SAf=5-j=d i ^/gg i >id j )f 

= d l r + A5r + [^{5g ir )d j \ r. 

The first term in the last line, being proportional to d l r, vanishes once a scalar product with n is taken. For the 
same reason, the only contribution that survives in the last term is obtained when the derivative di acts on f. The 
variation Sg^ is expressed in terms of Sr by taking the variation of g ik gk.j — 6*j so that finally: 

SH = h ■ A Sr — 2 (d l 8rd 3 r) (h ■ did r) . 

The variation of Ti c ,i = K i J d D SyfgH w.r.t. Sr involves the variation 5H and the variation of y/g, which is evaluated 
as before. The functional derivative of the first term is evaluated from the above equation with a partial integration: 
1 5H Cil 



= h\-HAf- d l fdiH + An + 2^-d t [y/g(h ■ d'djf^&r] \ 



V9 Sr 

On multiplying the above equation by h to find Q c ,\, the property n ■ dif = can be used again to show that the 
second term gives no contribution and the last term becomes 2b l jb J i . This is twice the trace of the square of the 
matrix of the coefficients of the second fundamental form. Finally we have to compute n- An. Using the facts that 
h -L d l h and that dih = —b^djf, we easily find: 



D 



h ■ A h = -d t h ■ d l h = b/ir (Dl) 



Collecting the various terms wc get the result displayed in Eq. (33) that is clearly fully R-invariant. 

The chain rule of differentiation, applied to Tt c ,p — K p J d D Syfg~H p , also gives Eq. (|35|), where the second term, as 
before comes from the variation of y/g while the others come from SH P = pH p ~ 1 SH. This also needs the straightfor- 



ward generalization of Eq. (Dl) to: 

h-A(Fh)= (a-^xAf 



for a generic R-invariant function F(s). 



APPENDIX E: ORIENTATIONAL ENERGY 



The variation of (|38|) for a change r — > r + Sr can be obtained as soon as we know S^fg and Sn. The former has 
already been obtained in appendix C and it is Sy/g = Jg gijdif '■ djSr whereas the latter is derived as follows. Since 
Sn _L h and 
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one has: 



= S(h ■ dif) = di? ■ 5n + h ■ di5f, 



Sh = —d l r(n ■ di5r). 



(El) 



The variation of Eq. (Bq) is then 



SH Z = I d D S ^g[-g^(d l r-d J Sr)x(n z ) 
, dx(n z 



dn. 



-{z ■ d l r) {h ■ di5r)] 



(E2) 



from which it follows that: 



i sn z 1 

= —On 



s[9 Sf Jg 3 

= X{n 
d\(n 



i d 2 x(n z ) dx(n z ) 

X{n z ) Ar — — n(z ■ d l r)d t n z niz ■ A rj 

dn z 



dn z 



dn 2 

-{d^r ■ djn z — (z ■ d->f)djn) 



(E3) 



The last two terms sum to zero since dih — —b^djf, and due to the symmetry of bij. Finally, using Eq. (A7) ( i.e., 
Ar = fiH) we get: 



1 sn z 

\/9 Sf 



nG z , 



with 



dx(n z )\ , d 2 X (n z ) , 
dni 



where z = z ■ f . 

Using Monge parametrization, 



r , n i dx{n z ) \ d 2 x(n z ) d^g 1 > 2 



dn 2 



with n z = g- 1 / 2 and H = VF^| 



(E4) 



(E5) 



(E6) 
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